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SUMMARY 

A predictive method for the titled flows based on the Prandtl energy method is 
developed and assessed by comparing predicted results with experimental results. 

For constant -density flows, both gross properties such as spreading rate and maximum 
turbulent kinetic energy and detailed properties such as mean shear stress distributions 
are shown to be well predicted. For variable -density flows, considerable attention is 
devoted to the inclusion in the analysis of the added effect of pressure fluctuations and 
of the variation in the several extant empirical parameters on the turbulent kinetic 
energy. It is found that a variation with Mach number of the characteristic Reynolds 
number for turbulent transport is needed to account for the observed decrease in 
spreading rate. The predictions which result from these considerations are compared 
with the limited experimental data presently available for the two crucial cases: com- 
pressible adiabatic mixing and low-speed isothermal mixing of two dissimilar gases. 

INTRODUCTION 

There is considerable activity presently underway by several groups on the 
description of turbulent shear flows by methods which contain more of the physics of 
turbulence than do the well-known methods based on mixing length and the usual eddy- 
viscosity models. Bradshaw et al., Rodi and Spalding, Hanjalic and Launder, Donaldson, 
and Wilcox and Alber (refs. 1 to 5) provide recent entry points in this literature. Most 
of this activity has concerned turbulent flows with uniform properties, presumably 
because these new methods should first be assessed for this simpler case and because 
the amount of data for the case of variable density is limited. Concern here is with the 
simplest turbulent shear flow, the two-dimensional mixing of two different streams, but 
under circumstances involving variable density, either because of the high speed of one 
stream or because of different compositions of the two streams. 

The two-dimensional mixing layer has the virtue that its description is given in 
terms of a similarity variable; thus, the numerical analysis associated with its study is 
modest. It has the further advantage that there are no solid walls present so that most 
effects of molecular transport are negligible and the modeling required to effect closure 


in the describing equations is simplified. A disadvantage is that its study in the labora- 
tory is relatively difficult when compared with other free mixing flows, for example, 
wakes and jets; therefore, experimental data suitable for purposes of comparison are 
limited, especially for the case of variable density. 

Previous theoretical work on two-dimensional, turbulent mixing for constant fluid 
properties is extensive; the well-known references 6 and 7 may be examined. For the 
case of variable density due to either compressibility effects of high speed or to hetero- 
geneous composition, the literature involving a serious effort to compare prediction and 
experimental data is sparse; most entries set up a model for the case of variable proper- 
ties, validate it against the data for the constant -density case, and then use it to predict 
the effects of variable density. Typical thereof is reference 8. 

This situation is perhaps due to the disarray in the experimental data presently 
existent in two crucial cases of two-dimensional mixing with variable density. The data 
have been recently reviewed by Birch and Eggers (paper no. 2 of this conference). With 
respect to the case of one high-speed stream mixing with a quiescent gas of the same 
composition under conditions such that the stagnation temperature is constant everywhere, 
that is, the case of so-called "compressible adiabatic flow," several sets of data indicate 
no effect of high speed (that is, of high Mach number) on the spreading angle. Other data 
indicate a significant increase in spreadii^ parameter, that is, a decrease in spreading 
angle and mixing rate, as the Mach number of the high-speed stream increases. 

The other crucial case of two-dimensional mixing with variable density involves the 
low-speed isothermal mixing of two gases of different molecular weights. Again, Birch 
et al. have pointed out the inconsistency of several sets of data related to the spreading 
parameter , in these cases to be considered a function of the velocity and density ratios 
of the two participating streams. 

Given this situation, any theoretical analysis such as the present one must be 
treated as a provisional one until at least the experimental evidence related to these two 
crucial cases cited is considered to be well-established. In the present work, the 
Brown and Roshko data (ref. 9) for heterogeneous, low-speed mixing and the data shown 
by Brown and Roshko and by Birch and Eggers indicating significant effect of Mach num- 
ber on the spreading parameter for compressible adiabatic flows have been accepted as 
correct. Accordingly, comparisons have been made of the present theoretical predic- 
tions against these data and the appropriate empirical constants have been adjusted to 
bring prediction and these experimental data into agreement insofar as possible. If future 
developments do not support these experimental data as correct, the present work should 
still provide a methodological framework of some value, but, of course, the adjustment of 
parameters would require some alteration. 
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A cursory examination of the new methods of analysis of turbulent shear flows indi- 
cates a wide variety of approaches, each of which can be extended to the case of variable 
density. One of the simpler ones, that is sometimes termed the "Prandtl energy method," 
is followed here. It is based on the idea of an eddy viscosity proportional to a length 
associated with the scale of the mixing layer and to the square root of the turbulent 
kinetic energy. It can be seen that the extension of this idea to flows with variable den- 
sity is ambiguous and must be guided by comparison with data. Comparison is made by 
relying heavily on the data associated with the gross behavior of the mixing layer , as 
reviewed and codified recently by Brown and Roshko (ref. 9) and Birch and Eggers (paper 
no. 2). 

The appropriate conservation equations are first developed and the terms as 
required for closure are modeled. One representation is then selected for the eddy vis- 
cosity for constant -density flows and the predictions based thereon are compared with 
various experimental results. The problem of incorporating the effects of variable den- 
sity in such a fashion as to achieve agreement with experimental data relative to the two 
crucial cases cited is then discussed. Finally, a comparison with experimental data is 
made. 


ANALYSIS 

The idealized flow shown schematically in figure 1 is considered. Two streams 
with different composition, velocity, and energy but with the same uniform pressure 
undergo turbulent mixing at the end of a splitter plate which is the origin of a x j ,X 2 
coordinate system. The symbols used in the analysis are defined in appendix A. 

In developing the equations for the description of this flow, several assumptions are 

employed. The effect of molecular transport is neglected except for certain dissipation 

terms. Constant pressure in the two external flows, flow similarity, and no chemical 

reaction are assumed. Furthermore, work is done in terms of mass -averaged quantities 

after Favre (ref. 10). For clarity a tilde will be used for a mass -averaged temporal 

mean; a double prime, for the fluctuating part left over; and a conventional bar, for a 

pu. 

regular temporal mean; thus, Uj(xpX 2 ,Xg,t) = + u^" = + u^" with customary nota- 

, _____ ^ 
tion. Note that u^" # 0 whereas puj" = 0. 

Basic Equations 

In terms of Cartesian tensor notation and in accord with these assumptions, it is 
readily established that the describing equations for the mean flow are 
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(i = 1, 2, . . .,N) 


* P“k"“i") = - 1: 
•^(pfSk^s * pVV) = 0 

iv 

^(puk?i + pUk”Yi") = 0 




( 1 ) 


a = 1, 2 N)j 


where is the mass fraction of species i and where hg is the stagnation enthalpy. 
In more detail 



Yihi + iukUk = h + lukUi 


k“k 


( 2 ) 


where h^ is the static enthalpy of species i. From equation (2) , 

pKg = ph + |(p'uj^ui^ + pV’V’) = + I ^(^k^k + 

where the mass -averaged turbulent kinetic energy is defined as 


q2 - pyv 


(3) 


Thus hg contains all forms of energy of present interest. 

For this particular closure scheme based on an eddy viscosity as discussed, a con- 
servation equation for q2 is needed; it is found to be 


0 PUj."u."u.’') + pu."u, " « -u, " - T., 

2 ax^V K ^ k 1 1 / ^ 1 k axjj k 9xi, ax,. 


(4) 


where is the viscous stress tensor. 


Closure Assumptions 

The closure of these equations as applicable to a thin shear layer wherein the 
boundary -layer approximations apply is next considered. First, in all the applications 
considered herein, it appears adequate to assume a single eddy transport coefficient 
essentially based on the mean velocity field. The generalization to a separate coefficient 
for each species, for the energy, and so forth is straightforward if an appropriate, sepa- 
rate lei^h scale for each coefficient is introduced. Additional comments on this will be 
made. Thus for the present a single eddy transport coefficient is introduced and 
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( 5 ) 



The remaining terms to be modeled arise in the equation for the turbulent kinetic 
energy. Dependence rests heavily on previous work related to the new phenomenology of 
turbulent shear flows; previous ideas are adapted to the a priori and formal introduction 
of an eddy viscosity. Thus for the dissipation term, 


8ui” 


axi, 



( 6 ) 


where is a constant for flows with constant density but is considered provisionally a 
function of an appropriate density ratio for flows with variable density. Equation (6) may 

be compared with the usual form ^ik(®“i"/^k) where L is a suitable length 

scale. Here dimensional arguments are used to let L «: e/q but a density ratio pjp\ 
raised to some power could be introduced without compromising dimensionality. 

It is customary in previous work devoted principally to constant -density flows to 
group the triple correlation terms and the pressure -velocity correlation together. What 
remains is a pressure -strain correlation which for constant -density flows is zero. Here 
more care must be taken; in appendix B the maintenance of the customary grouping and 
the addition of several new terms to account for the pressure -strain correlation prevail- 
ing in variable -density flows are heuristically justified. There results (see eq. (B9)) 




+ u. 






. (I . (,) 

where ^2 functions of a density ratio, otherwise constants, and <p 

is a thermodynamic parameter. (See appendix B.) This procedure should be compared 

0 / 0 \ ^ 

with the usual form for constant density, that is, — — (qL Here again L in terms 

9X2 V ®2/ 

of € has been eliminated and the terms appropriate for variable density have been 
added. 
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The mean pressure p is now considered; the Xg -momentum equation with the pres- 
sure in the external flows set to zero yields p = -pu2"u2". Thus in the Xj -momentum 
equation, 

^(p + pUi"U2’') = - PU2"U2") 

Although it is recognized that unless the turbulent kinetic energy is equally distributed 
amor^ the three velocity components, this term is nonzero; it does not appear to con- 
tribute significantly and will be dropped. 


Similarity Form 


The describing equations are now transformed to similarity form and the nondimen- 
sional variables are introduced. It is considered to be convenient to include a density 
distortion of the X2 -coordinate since the resulting equations are formally simpler and 
since the inverse transformation back to the physical variables is readily performed once 
a similarity solution is obtained. Thus, let 


Xi J( 


P A ^ 

dx« 

0 Pi 2 


and introduce 


f’ = 

U 
h 
h. 


H 

si 


Q2.ai 

U 2 




. IM.] 






( 8 ) 


(9) 


where a prime denotes differentiation with respect to r]. Thus, 


(*/')’ +«" = 0 

( 6 qH’)’ +fH’ =0 



The first two of these equations imply a Crocco relation 


Ho - 1 

H=i+-£ (1-f) 


432 


where Hg = (hg2/^sl)* equations (10) may be considered two equations for the 

two unknowns f(7j) and Q(7?), provided the density ratio and the eddy -viscosity 

parameter Cq are appropriately given. 

The boundary conditions are 

f’(oo) = l, f’(.=o) = l_y, f(0)=0 


( 12 ) 


Q(±oo) = 0 

Several remarks as to the boundary conditions on t{ri) are perhaps indicated. 


First, it is known from the work of Tii^ (ref. 11) that in applications of analyses of 
laminar or turbulent two-dimensional mixing of the sort considered here, the actual loca- 
tion of the dividing streamline f = 0 depends on details of the external flow, the pres- 
ence of walls, and so forth. Thus, the coordinates Xj and Xg must be considered to 
be boundary-layer coordinates along the actual dividing streamline whose location in 
physical coordinates may be found for a particular flow situation by application of 
reference 11. 


Second, in the presentation of experimental data relative to two-dimensional mixing 
it is customary to select as the origin of Xg the line along which f’ has a particular 
value, for example, f*(0) = 1 “ The difference between these two means for select- 
ing an origin is simply a translation in rj so that comparison between the solutions and 
such data is readily possible. 


Form for €q 

Since the eddy transport parameter €q clearly plays a central role in these equa- 
tions, it is appropriate to make some remarks about it. In terms of an academic inves- 
tigation, devoid of connection with the available data, a wide variety of relations for Cq 
can be assumed. Indeed, even if experimental results are conscientiously considered, 
there remain several such relations and in the course of this study several have been 
examined for flows with constant and variable density. Only those which involve coupling 
between the mean velocity and the turbulent kinetic energy, that is, those consistent with 
the spirit of the Prandtl energy method are reported. 


Comparisons of prediction and experiment lean heavily on the gross property 
spreading rate, as specified by Brown and Roshko (ref. 9). Define a parameter o such 


that 


where 


L32 

[(* 2 Ai), * (-^’‘1)2] 

Xg^Xj^j , i = 1, 2 are defined 


by the velocity ratios 
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(13) 


ui - (1 - y)u 

yU 



Uj - (1 - y)U 
yU 


(0. 1)^/2 


(i= 1) 


(i = 2)J 


Note that in defining a in terms of the mass-averaged velocity rather than in terms 
of the usual Uj^ which is implicitly assumed to be the mean velocity given by experi- 
mental data, the density -velocity correlations near the outer edges of the mixing layer are 
neglected since = Uj 

Note also that there are other possible definitions of quantities defining the spread- 
ing rate but equation (13) is convenient and according to Brown and Roshko (ref. 9) yields 
values of a "very close to those obtained from the more elaborate . . definitions. 
Accordingly, in considering some experimental data not treated by Brown and Roshko, 
the values of a given by the experimentalist have been taken and have been assumed to 
be equal to that which would be given by equations (13) if the detailed velocity profiles 
had been available. For the case y = 1, the most often quoted value a = Oq = 11.3 is 
due to Liepmann and Laufer (ref. 12). However, Birch and Eggers (paper no. 2) have 
pointed out that many experimentalists assume that this value is the value they would 
measure if they, in fact, did so, and has provided a more rational means for estimating 
the value of Oq peculiar to their setup. 


In the special case of constant density. 


(14) 


where A is a scale length on the order of the mixing-layer thickness. Now in the spirit 
of the Prandtl energy method it is assumed that at least for constant density 


^ = Rq = Constant 


(15) 


that is, there is a characteristic Reynolds number whose exact value depends on how A 
is defined. In cases requiring several eddy transport coefficients, it would be hoped that 
Rq would be constant for each value, provided consistent definitions of A are applied 
to each property being mixed. 

For the single eddy transport coefficient based on the velocity and for the two- 
dimensional mixing considered here, it is convenient to choose for A a thickness defined 
by 

A = [U - U(1 - y^ 


9ui 


(16) 
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This procedure is reminiscent of one means for defining the thickness of a shock wave 
and is readily generalized to other properties being mixed. Thus, the form of €q 
examined for flows with constant density becomes 


— Z Q 

■0 Rqf"(0)^ 


(17) 


The predictions of the analysis using equation (17) are compared with a variety of data 
for mixing having constant density and then €q is reconsidered for the case of variable 
density. 

In connection with the comparison of experiment and prediction, it is convenient to 
identify the similarity variable ^ ^ ~ ~ 


CONSTANT -DENSITY FLOWS 

For constant -density flows with equation (17) taken to relate €q to the dependent 
variables, a single parameter y specifies the flow situation but Rq, and ^2 


must be selected. The parameter (p=l,and a is immaterial for these flows. For 
the following consideration is made: if in the last of equations (10), specialized to 



This result can be identified with the relation between the mean shear stress and the tur- 
bulent kinetic energy widely used in the new methods of analysis of turbulent shear flows. 
(See, for example, ref. 1.) The generally accepted constant in this relation suggests 
= 0.024, the value used. 

Furthermore, for ^ 2 ^ which appears not to be critical, it is found on the basis of 
numerical experimentation that an appropriate value is 0.5. Finally, again on the basis 
of numerical experimentation, Rq = 22 is taken. 

Each of these three parameters has not been systematically varied to establish a 
set which is "optimum" in some sense but rather it has been determined on the basis of 
the comparisons of experiment and prediction that the cited set is adequate for many pur- 
poses. Also note that equations (10) have been solved by finite -difference methods and by 
iteration to handle their nonlinearity, the problem of the three-point boundary conditions, 
and the implicit appearance in the differential equations of f"(0), a quantity obtained from 
the solution. No features of the numerical analysis appear to warrant special mention; in 
fact, the computer program was designed primarily to provide means for altering readily 
the several parameters and the various flow situations to be studied and not to provide 
efficient computation of any one case. 
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Now the predictions can be compared with experimental data. Solutions have been 
computed for a variety of values of y and several parameters are shown in table I 
which are obtained from these solutions and which may have residual value. Note that 
Oq = 11.4 is predicted. 

In accord with previous remarks, figure 2 shows the prediction for Oq/o as 
well as the analytic, empirical approximation due to Sabin (ref. 13), CTq/ct s y/(2 - y). 
This approximation was shown by Birch et al. to represent well a wide variety of data 
if CTq is appropriately estimated for those cases in which it is not directly measured. 

It is seen that the predictions for are in good agreement with this empirical 

approximation and thus with experiment. 

Next, consider a second gross property. Yule (ref. 14) has recently developed an 
empirical equation for the maximum turbulent kinetic energy as a function of y. In 
terms of the variables used herein and with Oq = 11.4 for consistency, his equation 
becomes 

Q^(0) = 0.054(2 - y)y2 ^8) 

In figure 3 the predictions of this p^er are compared with those of equation (18). Again, 
very reasonable agreement is obtained. 

Because of the apparent high quality and completeness of the measured details of 
Spencer and Jones (ref. 15), they have been used to make comparisons with the predic- 
tions of this analysis. In figure 4 for y = 0,7, the measured and predicted mean shear 
stresses in the form -eQf" =* ^Uj are compared; in figure 5 a comparison is 

made for the turbulent kinetic energy in the form q 2 » (uj^ /u2). In both cases the 
agreement is good. It, of course, follows from the agreement in figure 4 that the mea- 
sured and predicted mean velocity profiles (uj/u) will also agree well. There seems 
little point in actually showing this comparison. 

It is also perhaps of interest to compare the distribution of eddy viscosity inferred 
from the measurements of mean shear and mean velocity with that computed. Although 
the experimental results may be considered subject to possible error, the results of 
carrying out the comparison are shown in figure 6 for Xj = 23.5 inches (59.69 cm) (dis- 
tance from the virtual origin), and for y = 0.7. The agreement must be considered highly 
satisfactory. 

Thus on the basis of these comparisons, it is concluded that the analysis for two- 
dimensional mixing flows with constant density gives satisfactory agreement over the full 
span of velocity ratios with gross properties such as spreading rate and with some of the 
available detailed properties. Further refinements of the constants /3j, /Sg, and Rq 
could presumably be made but to do so is not the main purpose of the present work. 
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VARIABLE -DENSITY FLOWS 


In the case of variable -density flows the constants j3j, niust be sup- 

plemented with the parameter a and the thermodynamic parameter <p. It is perhaps 
appropriate to make some general remarks about the approach used to incorporate the 
variable -density effects into the analysis. 

Crucial Cases and Alterations of Analysis 

As mentioned in "Introduction" there appear to be two crucial cases of variable- 
density, two-dimensional turbulent mixii^. In the compressible adiabatic mixing, y = 1 
and to a good approximation 

-f'2) <19) 

p 1 - m 

where m= (uy2hg^) = ‘ ^ ^ calorically 

perfect gas, and where is the Mach number in the high speed stream.* Equa- 

tion (19) gives the density ratio P 1 /P 2 = (1 “ 

The second crucial case is the low-speed isothermal mixing of two gases of differ- 
ent molecular weights, Wi and W 2 . In this case, 

-f') (20) 

p y 

where w s In this case the density ratio for y = 1 is Pi^P 2 = w. 

The experimental data which apply to the first case, compressible adiabatic mix- 
ing, and which appears to be most reliable (see p^er no. 2 by Birch and Eggers) indicate 
that a significant reduction in a occurs as Mj^ increases; in particular , a » 38 for 
m = 0.833 (Mj = 5). On the contrary, the most reliable data for the second case seem 
to be that of Brown and Roshko (ref. 9) which indicate no significant change in a with w 
over the range 0.143 ^ w i 7. The essential data for y = 1 are shown in figure 7 in 
terms of <Jq/cf plotted against Pi/P 2 * On the basis of a comparison of equations (19) 
and (20), one would not expect such a disparate difference in mixing behavior. 


* In this relation for m, y is the ratio of specific heats; there should be no con- 
fusion with our other use of y. In deriving equation (19), the contribution of turbulent 
kinetic energy to the stagnation enthalpy is essentially canceled by the variation of static 

pressure due to pu^^^u^. It is perhaps appropriate to observe here that m is the 

proper similarity parameter for compressible adiabatic mixing of calorically perfect 
gases so that high-speed mixing of gases other than air, for example, high-speed helium 
with quiescent helium, should be correlated in terms of m. 
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In the interest of exposition it is indicated now that with the parameters /3j, 
and Rq equal to the values used for constant -density flows and with a = 0.25 and cp 
appropriately selected, the analysis results in predictions in reasonable agreement with 
those of Brown and Roshko (ref. 9). However, it underpredicts very significantly the 
effect of highspeed on the spreading parameter. Therefore, accepting the experimentally 
observed difference in the two crucial cases to be correct, that is, anticipating that it will 
be substantiated by further experimental results, the implications thereof are now dis- 
cussed, that is, determination of the effect of Mach number through the parameter m. 
This apparent difference, accepted herein, in behavior of the two crucial cases rules out 
several possible means for incorporating the effects of variable density. Several means 
have been studied but have been abandoned, at least for the time being. 

The most obvious candidate to account for the difference in the two cases is the 
pressure fluctuations which are expected to become more significant with Mach number. 

It is shown in appendix B that in the variable -density case, the pressure correlations 
introduce effects into the equation for the turbulent kinetic energy through the diffusion 
term, the convection term, and through an added term directly associated with the density 
gradient. However, this equation, and thus presumably the turbulent kinetic energy, is 
dominated by the production and dissipation terms, neither of which are directly influ- 
enced by the pressure fluctuations. Thus it is found by numerical experimentation that 
unreasonable values for ^ 2 > treated as a function of m, must be assumed in order to 
achieve a significant alteration of a with m. It thus appears that at least up to a Mach 
number of 5, pressure effects cannot account for a significant compressibility effect on 
two-dimensional turbulent mixing. This result is in accord with the finding of Laufer 
(ref. 16) for turbulent boundary layers. 

K the changes in the empirical parameters with m that might be rationalized to 
account for the observed effect on a are considered, the parameter which accounts 
directly for dissipation and Rq appearing in €q remain. Again by numerical experi- 
mentation, it is found that unreasonably large changes in with m are required to 
bring prediction and experiment into agreement. Moreover, there appears to be no way 
to rationalize such changes; in this regard, note that such would imply a considerable 
deviation with m of the coefficient in the frequently assumed relation between mean 
shear and turbulent kinetic energy. Although it is not certain that some alteration of 
this coefficient is not in fact indicated, the alteration associated with the required changes 
in /3j is probably excessive. 

Thus, alterations of Rq with Mach number, that is, of the only parameter at our 
disposal, are considered. As is seen in more detail subsequently, rather modest changes 
in Rq with m can bring the predictions in agreement with the data for compressible 
adiabatic mixing. Consider here the possible physical explanation therefor, and keep in 
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mind a distinction between density and Mach number effects. With all the other param- 
eters in €q fixed, the value of Rq determines the rate of entrainment of external flow 
into the turbulent layer. This entrainment is known to be related to the detailed mechan- 
isms connected with the superlayer at the interface between the potential and turbulent 
flows. Because of detailed experiments involving conditioned sampling (see Kovasznay 
et al. (ref. 17) and Kaplan and Laufer (ref. 18)), a great deal is known about these 
mechanisms for low -speed turbulent boundary layers, for example, about the relative 
speed of the large turbulent eddies and the external flow. Although the corresponding 
data for the high-speed mixing layer do not exist and will be extremely difficult to obtain, 
it appears intuitively clear that when the relative speed of the large-scale turbulent 
eddies and the external flow approaches a significant fraction of the local speed of sound, 
major changes in the detailed entrainment mechanisms could occur. The observed 
decrease in spreading angle and mixing rate with free -stream Mach number, for example, 
as implied by figure 7, supports this view. These effects are concluded to be attributed 
to Rq in the present analysis. 

Accordingly, a strategy of seeking Rq = Rq(m) is adopted to bring the analysis 
into agreement with the case of compressible adiabatic mixing. Beyond the case of 
y = 1, that is, of mixing with a quiescent gas, it is expected that Rq will depend on y 
as well as on m when the second stream becomes supersonic. Resolution of this effect 
awaits further experimental data. 

Isothermal Binary Mixing 

Consider now the low -speed isothermal mixing of two dissimilar gases and treat 
w s (Wj^2) velocity ratio by means of y as parameters. The assumed values 

of w include helium-air cases with the molecular weight of air assumed to be 28. The 
value of (p is taken for simplicity to be three and represents its arithmetic mean for 
air and helium; these results are independent of any reasonable value for cp since it 
enters only in the convection term, that is, in one of the small terms, in the equation for 
the turbulent kinetic energy. Also ot = 0.25 is taken in all cases since it too does not 
enter importantly. Several of the combinations of w and y are chosen to permit com- 
parison with the experimental results of Brown and Roshko (ref. 9) for helium and air 
mixing. 

Table n gives the principal gross results from this series of calculations. Several 
points are indicated there. For a given velocity ratio, that is, y, the effect of the 
molecular weight ratio w on the spreading parameter a is small. This result is in 
accord with the main conclusion of Brown and Roshko (ref. 9). Moreover, simple calcu- 
lations using the results in table II show that for a given w, the effect of y on or is 
closely given by the equation of Sabin (ref. 13) developed for constant -density flows. This 
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prediction may be useful in determining Oq for binary mixing flows after the suggestion 
of Birch and Eggers (paper no. 2). 

The results with respect to the spreading parameter a for y = 1 are shown in 
figure 7 in terms of plotted against = w in this case. The comparisons 

with the results deduced by Brown and Roshko (ref. 9) are seen to be good, as suggested 
earlier; that is, when no changes in the effective parameters /3j, ^ 2 > Rq are made, 
the spreading parameter is insensitive to density variations. 

Attention is now turned to some of the detailed results of Brown and Roshko (ref. 9). 
In figure 8 a comparison is made of the predicted velocity and density profiles for one of 
their experiments corresponding to y = 0.622, w = 0.143, that is, helium at a higher 
velocity mixing with slower air. It is seen that the velocity profile is well predicted but 
that the density profile reflects the single eddy transport coefficient assumed a priori in 
the analysis, whereas the experimental results indicate a considerably larger coefficient 
for the species. Contrast this result with another case shown in figure 9, y = 0.622, 
w = 7, corresponding to faster moving air mixii^ with slower helium. Within the ability 
to read the plots in reference 9, it is not possible to distinguish the prediction from 
experiment so only the predicted results are shown. The reason for the difference in 
the predictability of the two cases is not clear. 

It is of interest to consider the predicted alteration in the distribution of some mean 
quantities with molecular weight or density ratio. Accordingly, the velocity profiles for 
the constant density case (w = 1) and the hydrogen-air cases w = 14, (1/14), each with 
zero velocity on one side of the mixing layer, are shown in figure 10. Although the 
spreading parameter is about the same in all three cases, it may be seen that the varia- 
tion in mean density results in altered velocity profiles. 

In figures 11 and 12 are shown the predicted distributions of production of turbulent 
kinetic energy in terms of eQf'’2(r/) and of dissipation in terms of 

plotted against ^ for these same three cases. There is a shift in the peaks of 

production and dissipation toward the low -density side of the dividing streamline and, as 
ejq)ected, the dissipation term is less than the production term in the central portion of the 
mixing layer. From the calculations it is found that the diffusion term provides most of 
the difference in the two main terms in the equation for the turbulent kinetic energy. It 
would be highly desirable if some measurements of the turbulent kinetic energy could be 
made in binary mixing flows to permit comparisons with predictions of the sort shown in 
figures 11 and 12 in order to establish whether the energy balances are as predicted or 
whether changes in the modeling are required. 

It is thus concluded that the analysis for constant -density flows without essential 
change ^the changes due to P^k'/^k) considered to be "unessential" here) largely 
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agrees with the results of Brown and Roshko (ref. 9) and shows little effect on the gross 
quantity, spreading parameter cr, although detailed distributions are as is to be expected 
altered by density variations. 


Compressible Adiabatic Mixing 


Attention is now turned to the case of compressible adiabatic mixing. Only y = 1, 
that is, mixing with a quiescent gas, is considered and thus there is a single parameter, 
that associated with the Mach number in the moving stream, fh. The previous values 
= 0.024, jSg = 0.5, and a = 0.25 are retained and (p = 7/2 is taken, corresponding 
to air. Again, the values of a and <p are unimportant. However, as discussed pre- 
viously, Rq is adjusted to bring into agreement prediction and experiment with respect 
to the spreading parameter o. 


As a result of preliminary numerical experimentation, it is found that a simple 
variation of Rq with m gives a variation of spreading parameter with 

^l/^2 “ reasonable agreement with the data shown in figure 7; thus, 

Rq = 22 + 5.2(1 - m)’^ (0 £ m S 0.7)' 


( 21 ) 


Rq = 39 


(0.7 i m < 1) 


This variation of Rq is used for the calculations discussed here. The implication of 
this result in terms of the previous discussion of the possible physical explanation of a 
change in Rq with Mach number appears to be that at a Mach number of about 3 the 
large-scale eddies in the mixing layer readjust so that further effects of Mach number 
on the spreading parameter are small and so that the characteristic Reynolds number Rq 
may be taken as a constant at higher Mach numbers. This explanation is, of course, to 
be considered a conjecture based on a careful consideration of the present analysis and 
the experimental data presently considered definitive. 

With the use of equations (21) a series of compressible adiabatic flows correspond- 
ing to a range of m have been computed; the principal gross results are given in 
table ni and the variation of the spreadii^ parameter in terms of with P\jp 2 

shown in figure 7. It is seen that the simple variation of Rq with m given by equa- 
tions (21) results in increases in cr with m in good agreement with the experimental 
data. 


In figures 13 to 15 are shown the same distributions of mean quantities as in the 
case of binary mixing in order to provide some indication of the effect of high speed on 
turbulent mixing. From figure 13 it can be seen that a significant change in the veloc- 
ity profile takes place when m increases from 0 to 0.667 (Mj from 0 to »3) but that 
little chaise seems to take place as m increases further despite the doubling of the 
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density ratio P 1 /P 2 ^ increases from 0.667 to 0.833. This behavior is also 

reflected in figures 14 and 15 which show that the distributions of production and dissipa- 
tion of turbulent kinetic energy for m = 0.667 and m = 0.833 are alike but are very 
different from the low -speed case (in = 0). 

It may also be noted from figures 14 and 15 that the peaks in production and dissipa- 
tion shift toward the low -density side of the dividing streamline as in the case of binary 
mixing. Also to be noted is the close balance between production and dissipation which 
indicates, as suggested earlier, that the other terms in the equation for conservation of 
turbulent kinetic energy are dominated by these two terms. 

The General Case 

On the basis of the results for the two crucial cases of turbulent mixing studied in 

some detail, the present analysis could be used with the parameters a 

fixed, with cp determined from thermodynamic considerations, and with Rq(5^ as 
given by equations (21) to make predictions of the properties of a variety of mixing flows 
involvii^ heterogeneous composition, high-speed effects, and nonadiabaticity. These 
have not been carried out in the present work because there appear to be no data with 
which to compare prediction and experiment. When such data are available, it would be 
of considerable interest to determine the extent of the agreement and/or disagreement. 

Of considerable value in assessing the a priori assumption of a single transport 
coefficient and the overall accuracy of the present analysis for these fundamental flows, 
that is, for two-dimensional turbulent mixing, would be cases of heterogeneous mixing 
under nonisothermal conditions, for example, heated helium mixing with relatively cold 
air. The greater the statistical detail constituting the data the more valuable such 
experiments would be. 


CONCLUDING REMARKS 

A simple turbulent flow, the two-dimensional mixing layer, involving significant 
variations in density due either to heterogeneity in composition or to compressibility 
effects associated with the high speed of one flow has been analyzed in some detail. One 
of the simpler of the new methods of analysis of turbulent shear flows, that usually 
termed the "Prandtl energy method," has been used. It involves relating the eddy trans- 
port coefficient to the turbulent kinetic energy and to an appropriate length scale of the 
large-scale turbulent motion. Thus, one equation in the second level in the hierarchy of 
describing equations, that describing the conservation of turbulent kinetic energy, must 
be added to the usual set of equations describing the variables of the mean flow. 
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For two-dimensional mixing layers with constant density it is found that the analy- 
sis provides satisfactory predictions for gross properties of the flow over the entire 
range of velocity ratios of the two streams and for detailed properties for the one case, 
that is, one velocity ratio, examined. 

With respect to flows involving variable density, there is discussed the existing 
experimental data for two crucial cases, the low -speed isothermal mixing of two dis- 
similar gases and the high-speed adiabatic mixing of air. It has been assumed provi- 
sionally that the results of Brown and Roshko, which show little alteration of the spread- 
ing rate with large density differences in the first case, and the results which show 
significant reduction in the spreading rate with Mach number in the second case are cor- 
rect, and that they will be confirmed by further experiments. With this assumption fixing 
the strategy to be followed, it is found that an alteration with Mach number of the empir- 
ical parameters entering the analysis is required to bring prediction and experiment into 
agreement. The analysis suggests that neither the modeling of pressure rate of strain 
nor alterations in two parameters, those relating to diffusion and dissipation in the equa- 
tion for the turbulent kinetic energy, can reasonably account for the observed Mach num- 
ber effect. However, it is shown that a modest change in the Reynolds number charac- 
terizing the turbulence reduces the spreading rate in accord with experiment. A sug- 
gestion is made of how this required change is related to the physics of the turbulence. 

It is emphasized that if further experiments do not support the assumed behavior of 
the two-dimensional mixing layer in the two crucial cases cited, the methodology of the 
present work may retain some value, but the strategy used to bring prediction and ejq)eri- 
ment into agreement would, of course, be altered. It is also emphasized that this study 
indicates the dangers of casually extending the new methods of analysis of turbulent shear 
flows to cases involving significant variations of density. 
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APPENDIX A 


Cp 

f 

H 

«2 

h 

hi 

hs 

L 

Mj 

m 

P 

Q 

q 

Ro 

Rq 

T 

t 


SYMBOLS 

coefficient of specific heat 
coefficient of specific heat of i species 
stream function (see eq. (9)) 
stagnation enthalpy ratio (see eq. (9)) 
stagnation enthalpy ratio, h^g^sl 
static enthalpy 

static enthalpy of i species 
stagnation enthalpy (see eq. (2)) 
length scale 

Mach number in stream 1 

Mach number parameter (see eq. (19)) 

static pressure 

nondimensionalized turbulent kinetic energy 

turbulent kinetic energy, — - — — 

P 

universal gas constant 
Reynolds number of turbulence 
static temperature 
time 
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u 


w 

Wi 

w 

Yi 

ot 

^ 1>^2 

y 


V 

A 

p 

I 

P 

Pi 


o 


x-component of velocity in stream 1 
Cartesian velocity components, i = 1, 2, 3 
molecular weight of mixture 
molecular weight of i species 
molecular weight ratio, 

Cartesian coordinates, i = 1, 2, 3 

mass fraction of i species 

fraction of turbulent kinetic energy due to U 2 

empirical constants 

parameter determining velocity of stream 2; ratio of coefficients of 
specific heat 

eddy viscosity 

nondimensionalized eddy viscosity 
transformed similarity variable (see eq. (8)) 
length scale 

kinematic viscosity coefficient 
similarity variable, 
density 

density in stream 1 
spreading parameter 
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<7q spreading parameter when uj = 0 

^ik viscous stress tensor 

(p thermodynamic parameter (see eq. (B7)) 
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APPENDIX B 


THE PRESSURE -STRAIN CORRELATION 


Here the modeling of the velocity -pressure gradient correlation appearing in equa- 
tion (4) is considered. The first step is to write 


ap _ 




auk” 


(Bl) 


For constant -density flows the second set of terms on the right-hand side is zero. Here 
the general case of a gas mixture described by a perfect gas law is considered. First 


P = 



(B2) 


where Rq is the universal gas constant and W is the mixture molecular weight. The 
conservation equations devoid of significant molecular effects are 





iK) * = f 


(i = l,2,3) 

> (B3) 


= “ <‘='-2 

The pressure relative to that in the two streams external to the mixing layer is mea- 
sured so that P^Xj, ±°o, Xg, t^ = 0. In addition, the stagnation enthalpy is defined as 

N 


hs . h + i ^ Yihi(T) + I 

i=l 

From mass conservation and the equation of state, it is 



easy to show that 


u, 

^ axj, dKy. 


(B4) 


(B5) 
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Next from equations (B3) and (B4), 


and 


where 


Thus 



P 






(B6) 


Except for certain dissipation terms which would enter if molecular effects were 
included in equations (B3), equation (B6) is exact. Drastic simplifications are now made. 
It hardly seems reasonable to add to the result of the crude modeling usually associated 
with the first set of terms on the right-hand side of equation (Bl) and complicated 
modeled terms arising from the second set. First note that Ro/CpW is to a good 
approximation constant because the model specific heat c„^Wi is roughly constant. 

Thus let 


Ox,. 


CpW 


- 1 


9p »i 8p \ 

u. + u, ’’ - p - — 


9x, 


Ox, 


Ox, 


Then by letting cp~^ = ^RQ^CpW^, it is found from equation (Bl) that 


Note now that the X 2 -momentum equation in the boundary -layer approximation and 
with the external pressure set equal to zero yields 

p = -pu2"u2" = -o!M^ (B8) 

where a is a parameter (taken to be constant) representing the fraction of the total tur- 
bulent kinetic energy in the U 2 -velocity fluctuations. Most low -speed data for turbulent 
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shear flows indicate a » 1/4 but there does not appear to be any data related thereto 
for variable -density flows. 

With equation (B8) , equation (B7) can be rewritten to yield 


V f -^K") + <1 - «’>“ a^(|) (B9) 


Several remarks about equation (B9) are appropriate. Equation (B9) reduces to 
that for constant density flows if (p = 1 so the usual model in this case is recovered. 

For compressible adiabatic flows, <p = y/{y - 1) where y is the ratio of specific heats. 
For binary isothermal flows, cp is really a function of concentration but an adequate 
approximation would appear to be an average value across the mixing layer. Finally, 
note that the inclusion of density effects in the velocity -pressure gradient correlation 
introduces two new terms (one directly related to the mean density gradient, and one to 
the convection of turbulent kinetic energy) and modifies the usual combination of the dif- 
fusion of turbulent kinetic energy and the pressure -velocity correlation. 
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TABLE I.- NUMERICAL RESULTS FOR CONSTANT DENSITY FLOWS 


y 

f (0) 

£"(0) 

Q 2(0) 

^0 

1 

0.578 

6.07 

6.42 X 10"2 

11.4 

0.8 

0.643 

7.18 

4.17 

16.8 

0.7 

0.680 

7.74 

3.22 

20.6 

0.6 

0.721 

8.31 

2.37 

25.8 

0.4 

0.808 

9.46 

1.06 

43.1 

0.2 

0.902 

10.6 

0.265 

97.0 


TABLE II.- NUMERICAL RESULTS FOR BINARY FLOWS 


y 

w 

f'(V = 0) 

f"(77=0) 

q 2( o ) 

Pi 

P (0) 

°^0 

1 

1 

0.578 

6.07 

6.42 xlO -2 

1 

11.4 


2 

0.624 

8.63 

6.92 

1.38 

11.2 


4 

0.664 

13.0 

6.88 

2.01 

11.6 


7 

0.694 

18.7 

6.82 

2.83 

12.1 


14 

0.729 

30.7 

6.70 

4.53 

13.2 


0.5 

0.544 

4.59 

6.94 

0.772 

10.9 


0.25 

0.505 

3.68 

6.92 

0.629 

11.1 


0.143 

0.487 

3.20 

6.87 

0.560 

11.1 


0.0714 

0.450 

2.79 

6.88 

0.489 

11.8 

0.8 

7 

0.744 

22.3 

4.40 

2.92 

17.2 

0.622 


0.794 

25.6 

2.67 

2.99 

24.8 

0.4 


0.861 

29.8 

1.11 

3.08 

42.5 

0.2 


0.928 

33.7 

0.279 

3.15 

91.0 

0.857 

0.143 

0.529 

3.45 

5.10 

0.529 

15.0 

0.622 


0.635 

3.93 

2.70 

0.497 

25.2 

0.4 


0.756 

4.42 

1.12 

0.478 

50.0 

0.2 


0.875 

4.93 

0.279 

0.466 

105 
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TABLE m.- NUMERICAL RESULTS FOR COMPRESSIBLE ADIABATIC MIXING 

[. = 1 ] 


m 

f’(7,= 0) 

nt?=o) 

q2(0) 

Pi 

P (0) 

*^0 

Rq 

0 

0.578 

6.07 

6.42 X 10"2 

1 

11.4 

22 

0.333 

0.610 

14.6 

4.12 

1.31 

19.8 

30 

0.5 

0.626 

21.6 

3.52 

1.61 

24.4 

32 

0.667 

0.649 

38.2 

2.77 

2.16 

31.6 

38 

0.75 

0.664 

54.3 

2.47 

2.68 

36.4 

39 

0.833 

0.684 

77.3 

2.43 

3.66 

38.4 

39 

0.875 

0.700 

99.3 

2.41 

4.57 

40.0 

39 
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Figure 1.- Schematic representation of the flow. 



i-y 


Figure 2.- Variation of spreading parameter with 
velocity ratio. Constant density flows. 





Figure 3.- Variation of turbulent kinetic energy on the dividing streamline 
with velocity ratio. Constant density flows. 



Figure 4.- Comparison of predicted and experimental mean shear stress. 

r = 0.7; constant density flows. 


455 





i 


Figure 5.- Comparison of predicted and measured turbulent kinetic energy. 

y = 0.7; constant density flows. 



c 


Figure 6.- Comparison of predicted and measured distributions of eddy viscosity, 
y = 0.7; constant density flows; = 23.5 inches (59.69 cm). 
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Figure 7.- Variation of spreading parameter in variable -density flows. 



{ 

Figure 8.- Comparison of predicted and experimental mean velocity and 
mean density profiles, y = 0.622; w = 0.143. 
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Figure 9.- Velocity and mean density profiles. 7 * 0.622; w = 7. 



-1.5 -1.0 -0.5 0 0.5 1.0 1.5 X 10 
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Figure 10.- Predicted velocity profiles for several density ratios. 

Binary mixing; r = 1 . 
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Figure 11.- Predicted distributions of production of turbulent kinetic energy 
for several density ratios. Binary mixing; y = 1. 



f 

Figure 12.- Predicted distributions of dissipation of turbulent kinetic energy 
for several density ratios. Binary mixing; y = 1. 


459 






Figure 13.- Predicted velocity profiles for several values of the Mach number 
parameter. Compressible adiabatic mixing; y = 1. 



Figure 14.- Predicted distributions of production of turbulent kinetic energy for several 
values of the Mach number parameter. Compressible adiabatic mixing; y = 1. 
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Figure 15.- Predicted distributions of dissipation of turbulent kinetic energy for several 
values of the Mach number parameter. Compressible adiabatic mixing; y = 1. 
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DISCUSSION 


J. Laufer ; I think Libby showed us very well the extent of our ignorance as far as density 
fluctuations in the turbulent flow are concerned. At this stage I would like to just ask one 
question - and maybe that question should be more properly addressed to Professor 
Roshko. Dr. Libby pointed out the difference between the conventional average and mass 
average velocity. I suspect that in many cases, the experimentalist measures a mass 
average quantity with a pitot tube. Professor Roshko, could you tell us, is it a mass 
average or a mean velocity that you are measuring? 

A. Roshko ; Well, that is a little hard to say. The velocities are measured by getting a 
pitot tube reading for pu^, but there is a little question of what you are reading with a 
pitot tube because you have the turbulence part of the term. When we do have an accurate 
measurement of the average density, we simply divide one by the other. What that aver- 
age means is hard to say. We estimate that it’s perhaps within 10 percent. 

J. Laufer ; I would like to suggest just on the basis of the simple Bernoulli's equation 
that one uses in obtaining the velocity that that value is, in fact, a mass average and not 
a conventional average. 

S. J . Kline ; On the same point, we have some experiments going on in which we are try- 
ing to vary the fluctuations and see what happens in a pitot under some controlled condi- 
tions. I don't want to go into details, but perhaps we can shed some light on your ques- 
tion. One more comment, we are getting ready to measure time averages by using laser 
Doppler techniques that we are developing, and then you will have a little more informa- 
tion on this point. 

P. A. Libby ; John, that’s exactly why I said, "if we take the experimenter’s word for it 
that he is giving us the usual average," then we are, in fact, comparing two slightly dif- 
ferent things. I understand, as we have talked before, about the fact that the experimenter 
may be wrong, but I do think that we have to take his word for it. 

J . Ito; We have been working with heterogenous shear flows for applications to the gas- 
gas auxiliary propulsion on the space shuttle system. We have obtained quite a bit of 
data in the past 6 to 8 months with cold flow and chemical reactions. One of the things 
that we have determined empirically was that we can best correlate mixing rates with 
pressure -gradient -type terms when you look at heterogenous flows of variable densities. 

In our case, we used gaseous hydrogen and gaseous nitrogen to simulate our oxygen, and 
we found that our best correlating parameters were really of the two differing 

jets, and the ratios of them. 
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